16.451 Lecture 21 Smallest nucleus: The deuteron Nov. 22, 2005 



Important because: 

• deuterium is the lightest nucleus and the only bound N-N state 

• testing ground for state-of-the art models of the N-N interaction. 



Electron scattering measurements: 
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Because the deuteron has spin 1, 
there are 3 form factors to describe 
elastic scattering: the "charge" ( & c ), 
"electric quadrupole" ( G q ) and "magnetic 
(G M ) form factors. (JLab data) 


Combined Data 
Intrinsic Shape of the Deuteron 
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Interesting feature: strong attractive np force, but 
a void in the center - the deuteron is hollow! ... Why? 




Interpretation of quantum numbers: 
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S — S n + s p 
J = S + L = 1 

7i — (+) (+) (— 1) L = + => L = 0, 2,4 


• Of the possible quantum numbers, L = 0 has the lowest energy, so we expect 
the ground state to be L = 0, 5-1 (the deuteron has no excited states!) 

• The nonzero electric quadrupole moment suggests an admixture of L = 2 

(more later!) 


introduce Spectroscopic Notation: 2S+1t 

with naming convention: L = 0 is an S-state, L = 1 is a P-state, L = 2: D-state, etc... 
-» the deuteron conf iguration is primarily 



Isospin and the N-N system: 



The total wavefunction for two identical Fermions has to be antisymmetric w.r.to 
particle exchange: 

^ total ~ space x $spin * X isospin 


Central force problem: 

V space (.r,e,t) = m Y lm ( 6, f) 

with symmetry (-1) L given by the spherical harmonic functions 


Spin and Isospin conf igurations: 

S = 0 and 1 = 0 are antisymmetric ; 
5-1 and 1=1 are symmetric 


3 c state can only 
c * be I = 0 ! 




Magnetic Moment: |u d = 0.857 jli n 


In general, the magnetic moment is a quantum-mechanical vector; it must be aligned 
along the "natural symmetry axis" of the system, given by the total angular momentum: 

jU ~ J 

But we don't know the direction of J , only its "length" and z-projection as expectation 
values: 

/j 2 \ = J (J+l); (J z ) = mj = ( -J.... + J ) 


in a magnetic field, the energy depends on mj via A E =-(jU»B ) = - g 7 B ju N 


Strategy: we will define the magnetic moment by its maximal projection on the 
z-axis, defined by the direction of the magnetic field, with mj = J 



Use expectation values of operators to calculate the result 



Calculation of ju : 




— 9 j J Mn 


Subtle point: we have to make two successive 
projections to evaluate the magnetic moment 
according to our definition, and the spin and 
orbital contributions enter with different weights. 


1. Project onto the direction of J: 


2. Project onto the z-axis with mj = J: 


A 

jli * J = ju cos a = 

//• J 

Jj(j+ 1) 
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g j J jU N = ju qos a qos ft 



jU • J 


1 
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Next, we need to figure out the operator for jd 



spin and orbital contributions: 
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We already know the intrinsic magnetic moments of the proton and neutron, so these 
must correspond to the spin contributions to the magnetic moment operator: 

/' P = + 2.79 /j n = g sp S jU N =^> g sp = +5.58 
Mn 1.91 /j n g sn S //jy 9 s,n 3.83 


For the orbital part, there is a contribution from the proton only, corresponding 
to a circulating current loop (semiclassical sketch, but the result is correct) 



For the deuteron, we want to use the magnetic moment operator: 



Details: 
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because m n = m p , and L is the total orbital angular momentum ! 

2. L = 0 in the "S-state" ( 3 Sj ) but we will consider also a contribution from 

the "D-state" ( 3 Dj) as an exercise 

3. The proton and neutron couple to S = 1, and the deuteron has J = 1 

— ~2 9 — 2~ { (^ s >P S P 9s, n ^ n 

Trick: use S n + S p = S and write the operator as: 

= [^ (9s,p + 9s,n) $ + y (0s,p “S?fPr S p _ 5„) + yi)/% 

- - 

40 * 

But the proton and neutron spins are aligned, and 
so the second term has to give zero! 






continued.... 


So, effectively we can write for the deuteron: 

P = \ (((0s,p + 9s,n) S+ L).j) /J n 

Trick for expectation values: 

J = L + S; ( j2 )= J ( J + !) 

(j 2 \ = ^(L + S).(L + S)\ = (i L 2 + S 2 + 2 L.s\ 
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(s.J) = (s.L + s.s) = ^J 2 - L 2 - S 2 ^j + (s 2 ^ = j (^J 2 -L 2 + S 2 ^ 
^L.J^ = (l.L + L.s'j = +L 2 -S 2 j 


M(%) = 9 s ,p+9s, n ) Mn = M P +M n = 0.880 Mi v 

M ( D\) = ^;( 8 ~(9 s,p + 9 s,n)) M n = 0.310 /j n 



Comparison to experiment: 


/u d = 0.857 /j n 


This is intermediate between the 
S-state and D-state values: 


ju = 0.880 jU N 

/d ( 3 D 1 ) = 0.310 /j n 
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Suppose the wave function of the deuteron is a linear combination of S and D states: 





with 




Then we can adjust the coeff icients to explain the magnetic moment: 

ju d = (1 -b 2 ) + b 2 // ( 3 D 1 ) 



b 2 = 0.04, or a 4% D-state admixture accounts for the magnetic moment ! 



Deuteron Summary so far 11 


Quantum numbers: (J 71 , 1) = (1\ 0) favor a 3 5 1 conf iguration (with S = 1, L = 0) as 

the lowest energy n-p bound state 


Magnetic moment: ju = 0.857 jli n is 2.6% smaller than for a pure 3 5 1 state and is 

consistent with a linear combination of L = 0 and L = 2 components: 




Vd) = a 

\) + b 

3 D ^ with a 2 + b 2 =1 and b 2 = 0.04 



If this is correct, we can draw two conclusions about the N-N force: 


deuteron: 


\ 

Also: the "hole" in the middle of the deuteron means that at very short distances 
the N-N potential must be repulsive - the neutron and proton do not overlap! 



1. The S = 1 configuration has lowest energy (i.e., observed deuteron 
quantum numbers) , so the N-N potential must be more attractive 
for total spin S = 1 than for S = 0. 

2. The lowest energy solution for a spherically symmetric potential is 
purely L = 0, and the L = 2 wave functions are orthogonal to L = 0 wave 
functions, so a physical deuteron state with mixed symmetry can only 
arise if the N-N potential is not exactly spherically symmetric! 




Deuteron Quadrupole moment etc... 12 


Formalism for electrostatic moments: 

{see also b.J. Griffiths, ' Introduction to E I ectody ramies' Ch. 3, and 16.369) 

• Recall that an electric charge distribution of arbitrary shape can be described by an 
inf inite series of multipole moments (spherical harmonic expansion), with terms of 
higher L ref lecting more complicated departures from spherical symmetry 

• The energy of such a distribution interacting with an external electric field E = - VV 
can be written as: 


dV 

E in, = V(0) q + — 

i e 2 v 
Pz „ „ 2 

Qzz **• 

oz 

o 4 3z 2 

0 


where the charge distribution is located at the origin and has a symmetry axis along z; 
q is the electric charge or ' monopole moment' 
p z is the electric dipole moment 
Q zz is the electric quadrupole moment 

and in general, higher order multipole moments couple to higher derivatives of the 
external potential at the origin. 



Formalism, continued.... (axially symmetric case) 13 


The multipole moments are defined in general by: 


E f = J p(r) E e d 3 r 


with corresponding multipole operators proportional to 
spherical harmonic functions of order b 


E, = c, r' Y lo (0) 


The f irst few multipole operators and associated moments are: 


order, i 

Moment 

Symbol 

operator, E i 

0 

monopole 

R 

1 =V4 xr°Y 00 (9) 

1 

dipole 

Pz 

z = r cos 0 = r 1 Y 1O (0) 

2 

— 

quadrupole 

Q zz 

r 2 ( 3 cos 2 <9 1) = ^ r 2 Y 20 (9) 



Connection to nuclei and the deuteron problem: 
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For a nuclear system, the multipole moments are expectation values of the multipole 
operators, e.g. the electric charge: 

z 

q = J p(r) d 3 r = J y/*(r) E 0 y/^r) d 3 r = + Ze 

i= 1 


The electric charge density is proportional to the probability density, i.e. 
the wave function squared, summed up for all the protons in the system !!! 



We already know the electric charge of the deuteron, but what about higher moments? 




continued... 
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Basic symmetry property: no quantum system can have an electric dipole moment if its 

wave function has definite parity, i.e. if I'Fl 2 is even. 

(Incidentally, the possibility of a nonzero electric dipole moment of the 
neutron is of great current interest: so far the measured upper limit 
is |p n | < 10 -25 e - cm ....) 

In fact, all odd multipole moments must vanish for the same reason 

First nontrivial case: Electric Quadrupole Moment Q zz (F&Hsec. 18 . 1 ) 



z 




example: uniform density object with ellipsoidal shape 


Q > 0 if b > a: "prolate" 
Q < 0 if b < a: "oblate" 



How to relate this to a nucleus, e.g. the deuteron? 
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Assume that the charge distribution is an ellipsoid of revolution with symmetry axis 
along the total angular momentum vector J: 



As for the magnetic dipole moment, 
we specify the intrinsic electric 
quadrupole moment as the expectation 
value when J is maximally aligned with z: 


/ rv \ 


Qint = (^2 J 

i — mm 

3 

II 


n m J J 

"Quantum geometry": cos p - yyj - ^ 


If an electric field gradient is applied along the z-axis as shown, the observable 
energy will shift by an amount corresponding to the intrinsic quadrupole moment 
transformed to a coordinate system rotated through angle p to align with the z-axis 



Result: 
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Let Q| ab be the electric quadrupole moment we measure for the ellipsoidal charge 
distribution with mj = J: 



How do we apply this to anything? 

1. A spherically symmetric state (L = 0) has Q jnt = 0 (e.g. deuteron S-state) 

2. Even a distorted state with J = i will not have an observable quadrupole moment 

3. J = 1, L = 2 is the smallest value of total angular momentum for which we can observe 
a nonzero quadrupole moment 

-> these are the quantum numbers for the deuteron D- state! 


Deuteron intrinsic quadrupole moment: 


recall our model for the deuteron wave function: 




a 
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result for the quadrupole moment: 


Q,„, = — \a'b\ ( r 2 ) $D 


note cancellation here 


/ 


10 


1 b 2 (r 2 ) 

/ DD 
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= + 0.00286 ± 0.00003 bn 


A good model of the N-N interaction can fit both the magnetic moment and the 
quadrupole moment of the deuteron with the same values of a and b! -> 



Basic features of the N-N potential, via the deuteron, etc.: 
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1. Independent of the value of L, the state with intrinsic spins coupled to S = 1 
has lower energy 

-> this implies a term proportional to: 

' -1/4, S= 1 



+ 3/4, S=0 

V 

2. The deuteron quadrupole moment implies a non-central component, i.e. the potential 
is not spherically symmetric. Since the symmetry axis for Q is along J, Q > 0 
means that the matter distribution is stretched out along the J - axis: 




N-N interaction continued.... 
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3. There is also a spin-orbit term, as deduced from N-N scattering experiments with 
a polarized beam: 

V s _o - (L.S) 


(This plays a very important role also in determining the correct order of energy 
levels in nuclear spectra - more later!) 


4. Finally, all contributions to the N-N interaction are based on a microscopic meson 
exchange mechanism, as explored in assignment 5: 





Where M is a n, p, co ... meson, etc. 

and each term has a spatial dependence 
of the form: 

—m r 

e 

V(r) = g x(spin function ) 

r 



State of the art N-N interaction model: (not to scare anybody...) 21 


PHYSICS REPORTS (Review Section of Physics Letters) 149, No. 1 (1987) 1-89. North-Holland, Amsterdam 


(89 page exposition of one of only ~3 state-of-the-art models of the N-N 
interaction worldwide - constantly refined and updated since first release.) 


THE BONN MESON-EXCHANGE MODEL FOR THE 
NUCLEON-NUCLEON INTERACTION* 

R. MACHLEIDT* 

Los Alamos National Laboratory, MS H850, Los Alamos, NM 87545, U.S.A. 

and 

Department of Physics, University of California, Los Angeles, CA 90024, U.S.A. 

and 

K. HOLINDE ++ and Ch. ELSTER 

Institut fur Theoretische Kernphysik der Universitdt Bonn, Nussallee 14-16, D-5300 Bonn, West Germany 



Brace yourself - it takes 2 pages just to write all the terms down!!! 
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pseudoscalar mesons: 


V ps( m ps’ r ) = To 


1 £ps 


12 4t r mps 


m 


ps 


m 


f("V)o , i:®'2 + z (»v)s 


ps / 12 




(tensor operator) 


(F.6) 


scalar mesons: 


(they use "a" for spin) 


£ s 


V,(m s ,r) = - ^m s 


l-ip 

4 V m , 


2 -| 


1 r-r2i 


+ Y(«/)r 2 ] + \Z,{m s r)L • S f ; 


vector mesons: 


/ (F.7) 

(spin-orbit interaction) 


KK> 0 = 7 Z 


h. J 

4tt my l 


1+ T ( — 


2 ~\ 


Y(m v r) - — 2 [V 2 Y(m v r) + Y(m v r)V 2 ] 


+ - 


1 l m. 


+ -T 


6 \ m 

1 Sv/v 


) Y(m v r)or • o- 2 - f Z^m^L • S - i Z(m v r)S 12 


2 4tt 


m v {(m v /m) 2 Y(m v r) + |(m v /m) 2 Y(m v r)or • <r 2 


-4Zj(m v r)L • S - 5 Z(m v r),S 12 } 

f 2 

+ 4^ w v {K ^ v /m ) 2 y ( m v r )»'i ‘ ^2 “ AZ(m v r)S 12 } , 



more of this... 
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with 


Y(x) = &~ x /x , 


Z( x) = (mjm) 2 ( 1 + 3/x + 3 lx 2 )Y(x) , 



"Yukawa functions" 
and derivatives... 


Z 1 (x) = 



2 

( 1 lx + 1 lx 1 )Y{x) , 



(F.8) 


and 



We use units such that h = c = 1 (he = 197.3286 MeV fm). The use of the form factor, eq. (3.3), at each 
vertex (with n t( = 1) leads to the following extended expressions: 


K.W = K(m a , r) - 


Ala ~ m l 

A la ~ Ala 


K(4, i,r) + 


A l , i ~ 

Ala ~ <i 


V a (\a > . 


(F.9) 


where A. , = A, + s, = A, — s, e = 10 MeV is an appropriate ch oice 


The full NN potential is the sum of the contributions from six mesons: 
Hr) = 1 V„(r) 

a = / iT,p,T),a),8,a 


\ 



Parameters and Predictions: 
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R. Machleidt et ai, The Bonn meson-exchange model for the nucleon-nucleon interaction 


Table 14 


potential (OBEPR) 

deuteron properties: 

[fjg J 

m a (MeV) A. (GeV) 

LEP 

Theory 

TT 14.9 

138.03 

1.3 

e d (MeV) , 

2.2246 




Pn (%) 

4.81 

p 0.95; [6.1] 

769 

1.3 

2d (fm 2 ) 

0.274 




Md (f%) 

0.8524 

n 3 

548.8 

1.5 

-4 S (fm I;I ) 

0.8860 




D/S 

0.0260 

o) 20; [0.0] 

782.6 

1.5 

G (fm) 

1.9691 




a s (fm) 

-23.751 

5 2.6713 

983 

2.0 

r._ (fm) 

2.662 




a t (fm) 

5.423 

a 7.7823 a 

550 8 

2.0 

r, (fm) 

1.759 

Impressively good agreement for ~ 

10,000 experimental 


low energy scattering 

data points in assorted n- 

-p and p-p scattering 


parameters, etc: 



What does the NN potential look like? 25 


It looks different in different spectroscopic states of the 2N system! 



Only the deuteron is bound! Its quantum numbers have the deepest potential well. 




